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We develop a theory of superconductivity in ultrasmall (nm-scale) metallic grains having a discrete 
electronic eigenspectrum with mean level spacing d ~ A (bulk gap). The theory is based on 
calculating the eigenspectrum using a generalized BCS variational approach, whose (qualitative) 
applicability has been extensively demonstrated in studies of pairing correlations in nuclear physics. 
We discuss how conventional mean field theory breaks down with decreasing sample size, how 
the so-called blocking effect (the blocking of pair-scattering by unpaired electrons) weakens pairing 
correlations in states with non-zero total spin (thus generalizing a parity effect discussed previously) , 
and how this affects the discrete eigenspectrum's behavior in a magnetic field, which favors non-zero 
total spin. In ultrasmall grains, spin magnetism dominates orbital magnetism, just as in thin films 
in a parallel field; but whereas in the latter the magnetic-field induced transition to a normal state 
is known to be first-order, we show that in ultrasmall grains it is softened by finite size effects. Our 
calculations qualitatively reproduce the magnetic-field dependent tunneling spectra for individual 
aluminum grains measured recently by Ralph, Black and Tinkham [Phys. Rev. Lett. 78, 4087 (1997)]. 
We argue that previously-discussed parity effects for the odd-even ground state energy difference 
are presently not observable for experimental reasons, and propose an analogous parity effect for 
the pair-breaking energy that should be observable provided that the grain size can be controlled 
sufficiently well. Finally, experimental evidence is pointed out that the dominant role played by 
time-reversed pairs of states, well-established in bulk and in dirty superconductors, persists also in 
ultrasmall grains. 



PACS numbers: 74.20.Fg, 74.25.Ha, 74.80.Fp 
I. INTRODUCTION 

What happens to superconductivity when the sample 
is made very, very small? Andersorjj addressed this ques- 
tion already in 1959: he argued that if the sample is so 
small that its electronic eigenspectrum becomes discrete, 
with a mean level spacing d = 1/N(ef) ~ 1/Vol, "su- 
perconductivity would no longer be possible" when d be- 
comes larger than the bulk gap A. Hcuristically, this is 
obvious (see Fig. [l] below): A/d is the number of free- 
electron states that pair-correlate (those with energies 
within A of ep), i.e. the "number of Cooper pairs" in 
the system; when this becomes ~1, it clearly no longer 
makes sense to call the, system "superconducting" . 

Giaver and ZellerBtil were among the first to probe 
Anderson's criterion experimentally: studying tunneling 
through granular thin films containing electrically insu- 
lated Sn grains, they demonstrated the existence of an en- 
ergy gap for grain sizes right down to the critical size es- 
timated by Anderson (radii of 25 A in this case), but were 
unable to prove that smaller particles are always normaL 
Their concluding comments are remarkably perspicuous:™ 
"There can be no doubt, however, that in this size region 
the bulk theory of superconductivity loses its meaning. 
As a matter of fact, perhaps we should not even regard 
the particles as metallic because the energy-level spacing 
is large compared to kT and because there are very few 
electrons at the Fermi surface. The question of the lower 
size limit for superconductivity is, therefore, strongly cor- 
related with the definition of superconductivity itself." 

These remarks indicate succinctly why the study of su- 



perconductivity near its lower size limit is of fundamental 
interest: the conventional bulk BCS approach is not di- 
rectly applicable, and some basic elements of the theory 
need to be rethought, with the role of level discreteness 
demanding special attention. 

Eirst steps in this directionjvere taken by Strongin et 
al.u and by Muhlschlegel et alB, who calculated the ther- 
modynamic properties of small superconducting grains. 
However, since experiments at the time were limited to 
studying ensembles of small grains (e.g. granular films), 
there was no experimental incentive to develop a more de- 
tailed theory for an individual ultrasmall superconduct- 
ing grain, whose eigenspectrum, for example, would be 
expected to reveal very directly the interplay between 
level discreteness and pairing correlations. 

This changed dramatically in 1995, when Ralph, Black 
and Tinkham (RBT)Q succeeded in constructing a single- 
electron transistor (SET) whose island was an ultrasmall 
metallic grain: by studying the tunneling current through 
the device, they achieved the first measurement of the 
discrete eigenspectrum of a single grain. This gtLabled 
them to probe the effects of spin-orbit scattering. I a non- 
equilibrium excitations^ and superconductivity,Lra which 
manifests itself through the presence (absence) of a sub- 
stantial spectral gap in grains with an even (odd) number 
of electrons. 

RBT's work stimulated several theoretical investiga=. 
tions. Besides discussing non-equilibrium— effects JiJ'til 
these focused mainly on superconductivity,E2Tt£l and re- 
vealed that the breakdown of pairing correlations with 
decreasing grain size predicted by Anderson harbors 
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some-surprises when scrutinized in more detail: von Delft 
et al£3 showed that this breakdown is affected by the par- 
ity (p) of the number of electrons-ma-the grain: using 
parity-projected mean-field theorytZllia and variational 
methods and assuming uniformly spaced electron lev- 
els, they solved the parity-dependent gap equation for 
the even or odd ground state pairing parameters A e or 
A Q as function of d (using methods adapted from Stron- 
gin et al.u), and found that A (d) < A e (d), i.e. ground 
state pairing correlations break down sooner with in- 
creasing d in an odd than an even grain (the difference 
becoming significant fop-|C? ~ A). This is due to the 
so-called blocking effectxB the odd grain always has one 
unpaired electron, which blocks pair-scattering of other 
pairs and thereby-iweakens pairing correlations. Smith 
and Ambegaokart2l showed that this parity effect holds 
also for a random distribution of level spacings (as-also 
anticipated by BlantercJ) , and Matveev and Larkinli-a in- 
vestigated a parity effect occuring in the limit d 3> A. 

The A c < A e parity effect-has an obvious generaliza- 
tion, studied by Braun et alrB using a-generalized BCS 
variational approach due to Soloviev:£3 any state with 
non-zero spin s (not just the odd ground state) experi- 
ences a significant reduction in pairing correlations, since 
at least 2s electrons are unpaired, leading to an enhanced 
blocking effect (A s < A s / if s > s'). The latter's con- 
sequences can be observed in the magnetic-field depen- 
dence of SET tunneling spectra, since a magnetic field fa- 
vors states with non-zero spin and consequent enhanced 
blocking effect. In ultrasmall grains, spin magnetism 
dominates orJaital magnetism, just as in thin films in a 
parallel field^il but whereas in the latter the magnetic- 
field induced transition to a normal state is known to be 
first-order, Braun et al. showed that in ultrasmall grains 
the transition is softened due to finite size effects. More- 
over, they argued that some of RBT's grains fall in a 
region of "minimal superconductivity" , in which pairing 
correlations measurably exist at H = 0, but are so weak 
that they may be destroyed by the breaking of a single 
pair (since the number of electron pairs that take part in 
the formation of a correlated state becomes of order one 
for d ~ A). 

In the present paper we elaborate the methods used 
and results found by Braun et al. in Ref. [l^ and present a 
detailed theory of superconductivity in ultrasmall grains. 
Our discussion can be divided into two parts: in the first 
(sections || and III), we consider an isolated ultrasmall 
grain and (a) define when and in what sense it can be 
called "superconducting" ; (b) use a generalized BCS vari- 
ational approach to calculate the eigenenergies of various 
variational eigenstates of general spin |s), which illus- 
trates the break-down of mean-field theory; and (c) dis- 
cuss how an increasing magnetic field induces a transi- 
tion to a normal paramagnetic state. In the second part 
(section IV), we consider the grain coupled to leads as 
in RBT's SET experiments and discuss observable quan- 
tities: (a) We calculate theoretical tunneling spectra of 



the RBT type, finding qualitative agreement with RBT's 
measurements; (b) point out that the above-mentioned 
ground state energy parity effect can presently not be ob- 
served, and propose an analogous pair-breaking energy 
parity effect that should be observable in experiments of 
the present kind; and (c) explain how RBT's experiments 
give direct evidence for the dominance of time-reversed 
pairing, at least for small fields (implying that the suf- 
ficiency of using only a reduced BCS-Hamiltonian, well- 
established for bulk systems and dirty superconductors, 
holds for ultrasmall grains, too). 



II. PAIRING CORRELATIONS AT FIXED 
PARTICLE NUMBER 

The discrete energies measured in RBT's experiments 
essentially correspond to the eigenspectrum of a grain 
with fixed electr on nu mber N (for reasons explained in 
detail in section IV A ). In this and the next section, we 
therefore consider an ultrasmall grain completely isolated 
from the rest of the world, e.g. by infinitely thick oxide 
barriers. 

When considering a truly isolated superconductor (an- 
other example would be a superconductor levitating in 
a magnetic field due to the Meissner effect) one needs 
to address the question: How is one to incorporate the 
fixed-iV condition into BCS theory, and how important 
is it to do so? Although this issue is well understood 
and was discussed at length in the early days of BCS 
theory, in particular in its application to pairing corre- 
lations in nucleiE3 p ' 439 , for pedagogical reasons the ar- 
guments are worth recapitulating in the present context. 
We shall first recall that the notion of pair-mixingO that 
lies at the heart of BCS theory is by no means inherently 
grand-canonical and can easily be formulated in canon- 
ical language, then summarize what has been learned 
in nuclear physics about fixed-iV projection techniques, 
and finally conclude that for present purposes, standard 
grand-canonical BCS theory should be sufficient. Read- 
ers familiar with the relevant arguments may want to 
skip this section. 



A. Canonical Description of Pair-Mixing 

Conventional BCS theory gives a grand-canonical de- 
scription of the pairing correlations induced by the pres- 
ence of an attractive pairing interaction such as the re- 
duced BCS interaction 
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(with V > 0). (1) 



(The Cj± are electron destruction operators for the single- 
particle states \j, ±), taken to be time-reversed copies of 
each other, with energies Sj±.) The theory employs a 
grand-canonical ensemble, formulated on a Fock space of 
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states in which the total particle number TV is not fixed, 
as illustrated by BCS's variational ground state Ansatz 

\BCS) = Y[( Uj + vjc) + c\_) |Vac) (u 2 + v 2 = 1). (2) 



This is not an eigenstate of the number operator N — 



Yljr Cjo-Cja and its particle number is fixed only on the 
average by the condition (BCS\N\BCS) = N, which de- 
termines the grand-canonical chemical potential fi. Like- 
wise, the commonly used definition 
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(3) 



for the superconducting order parameter only makes 
sense in a grand-canonical ensemble, since it would triv- 
ially give zero when evaluated in a canonical ensemble, 
formulated on a strictly fixed- N Hilbert space of states. 

A theory of strictly fixed- N superconductivity must 
therefore entail modifications of conventional BCS the- 
ory. In particular, a construction different from Abcs 
is needed for the order parameter, which we shall hence- 
forth call "pairing parameter" , since "order parameter" 
carries the connotation of a phase transition, which would 
require the thermodynamic limit N — > oo. The pair- 
ing parameter should capture in a canonical framework 
BCS's essential insight about the nature of the super- 
conducting ground state: an attractive pairing interac- 
tion such as H re d will induce pairing correlations in the 
ground state that involve pair-mixing across Ef (see also 
Ref. [l2"| ), i.e. a non-zero amplitude to find a pair of time- 
reversed states occupied above ef or empty below ef- 
BCS chose to express this insight through the Ansatz 
(||), which allows Vj 7^ for Ej > £p and Uj =/= for 
£j < Ef- It should be appreciated, however (and is made 
clear on p. 1180 of their original papera), that they chose 
a grand- canonical construction purely for calculational 
convenience (the trick of using commuting products in 
(^J) makes it brilliantly easy to determine the variational 
parameters Uj, Vj), and proposed themselves to use its 
projection to fixed N, \BCS)n, as the actual ground 
state. 

Since [H re d, N] = 0, one would expect that the 
essence of BCS theory, namely the presence of pair- 
mixing and the reason why it occurs, can also be for- 
mulated in a canonically meaningful way. Indeed, this 
is easy: pair-mixing is present if the amplitude Vj = 
(c < j + c^_Cj_Cj + ) 1 / 2 to find a pair of states occupied is 
non-zero also for £j > £p, and the amplitude Uj = 
{cj-Cj+c^c^) 1 / 2 to find a pair of states empty is non- 
zero also for Ej < ef (the bars indicate that the Uj 
and Vj defined here differ in general from the Uj and 
Vj used by BCS; note, though, that the former reduce to 
the latter if evaluated using \BCS)). The intuitive rea- 
son why H rec i induces pair-mixing in the exact ground 
states \G) despite the kinetic energy cost incurred by 



shifting pairing amplitude from below to above Ef, is 
that this frees up phase space for pair-scattering, thus 
lowering the ground state expectation value of H re d'- 
in (G\H re d\G), the jf term can be non-zero only if 
\G) ^ 0, implying {v f ) G ^ and 



j+ c ]- c j'- c j'+ 



both c 

(uj)c 7^ 0, and also (Gr|ct + ct_c 7 -<_c J */ + ^ 0, implying 
{vj)c 7^ and (uji)a ^ 0. By pair-mixing, the system 
can arrange for a significant number of states to simul- 
taneously have both (vj)a 7^ an d {uj)G 7^ 0; this turns 
out to lower the ground state energy sufficiently through 
(G\H re d\G) that the kinetic energy cost of pair-mixing 
is more than compensated. Furthermore, an excitation 
that disrupts pairing correlations in the ground state by 
"breaking up a pair" will cost a finite amount of energy 
by blocking pair-scattering involving that pair. For ex- 
ample, the energy cost of having definitely occupied 
(v,j = 0) and \j— ) definitely empty (vj = 0) is 



(G|^ Cj W|G})+y(G| C t^, > r r ,_r 



\G), 



in which the restricted sum reflects the blocking of scat- 
tering involving the j-th pair. When evaluated using 
l-BCS*), this quantity reduces to Ej(l—2vJ)+UjVjABcs — 
[e| + Ages] 1 / 2 , which is the well-known quasi-particle 

energy of the state ^ + \BCS). 

The above simple arguments illustrate that there is 
nothing inherently grand-canonical about pair-mixing. 
Indeed, at least two natural ways suggest themselves to 
measure its strength in a canonically meaningful way, 



using for instance the pairing parameter A = 
proposed in Ref.ta, or one proposed by Ralphta: 
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(4) 



Both A and A' were constructed such that they reduce, 
as is desirable, to the same result as Abcs when each 
is evaluated using \BCS) with real coefficients Uj,Vj, 
namely Vj^j u j v j- A n appealing feature of A' is that by 
subtracting out (cj + Cj + ) (cj_Cj-) , it transparently em- 
phasizes the pairing nature of superconducting correla- 
tions, i.e. the fact that if \j+) is empty (or filled), so 
is A' will be very small if the occupation of \j+) 

is uncorrelated with that of as it is in a normal 

Fermi liquid. The overall behavior (as function of energy 
Ej) of the summands in both A and A' will be similar 
to that of UjVj (though not identical to UjVj or to each 
other; a quantitative evaluation of the differences, which 
increase with increasing d/A, requires an honest canoni- 
cal calculationEj) . UjVj is shown in Fig. 0(a), which illus- 
trates that pair-mixing correlations are strongest within 
a region of width Abcs- In this paper, we shall call 
a system "pair-correlated" if A' is a significant fraction 
of its bulk value A (say, somewhat arbitrarily, at least 
25%), and regard this as being synonymous with "super- 
conducting" . 
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B. On the breaking of Gauge Symmetry 

In some discussions of conventional BCS theory the 
defining feature of superconductivity is taken to be the 
breaking of gauge symmetry by the order parameter. 
This concept is illustrated by the BCS order parameter 
Abcs °f Eq. (||): if non-zero, it has a definite phase and 
is not gauge-invariant (under Cj a — > e l ^Cj ai it changes to 
e^Ascs)- Note, though, that this point of view cannot 
be carried over to fixed- TV systems. Firstly, these trivially 
have A bcs — 0, and secondly and more fundamentally, 
the breaking of gauge symmetry necessarily presupposes 
a grand-canonical ensemble: since phase and particle 
number are quantum-mechanically conjugate variables, 
formal considerations dictate that the order parameter 
acquire a definite phase only if the particle number is 
allowed to fluctuate, i.e. in a grand-canonical ensemble. 

Of course, in certain experimental situations where iV 
manifestly does fluctuate, such as the celebrated Joseph- 
son effect of two superconductors connected by a tun- 
nel junction, their order parameters do acquire definite 
phases, and their phase difference is a measurable quan- 
tity. However, for a truly isolated superconductor with 
fixed TV the "phase of the order parameter" is not ob- 
servable, and the concept of gauge symmetry breaking 
through an order parameter with a definite phase ceases 
to be useful. Indeed, the canonically meaningful pair- 
ing parameters A and A' defined above are manifestly 
gauge- invariant . 



C. Fixed-TV Projections 

It is easy to construct a variational ground state ex- 
hibiting pair-mixing and having definite particle number, 
by simply projecting \BCS) to fixed TV, as suggested by 
BCStE This can be achieved by the projection integral 



\BCS) N = f t 
Jo 



-i<pN 



IK* 



v 3 c\ + c\_) I Vac), (5) 



whose randomization of the phases of the Vj's illustrates, 
incidentally, why gauge invariance is not broken at fixed 
TV. 

This and related fixed-TV projections were studied in 
great detail in nuclear physics, with the aim of variation- 
ally calculating nuclear excitation spectra for finite nu- 
clei (TV < 240) exhibiting pairing correlations (Ring and 
Schuck provide an excellent review of the extensive liter- 
ature, see chapter 11 of Ref. |2^; a recent reference is 26). 
The simplest approach is called "projection after varia- 
tion": the unprojected expectation value (BCS\H\BCS) 
is minimized with respect to the variational parame- 
ters {vj}, which thus have their standard BCS values 
Vj = Ejj (ej+A^cg) 1 / 2 }, but then these are inserted 
into \BCS)n and expectation values evaluated with the 
latter instead of \BCS). This elimination of "wrong- TV" 



states after variation turns out to lower the ground state 
energy relative to the unprojected case (by a few percent 
in nuclei) and thus improves the trial wave-function. Fur- 
ther improvements are possible using the more sophisti- 
cated "projection before variation" strategy, where the 
projected expectation value n(BCS\H\BCS)n is mini- 
mized with respect to the {vj}. However, these then no 
longer have the simple BCS form, but instead are deter- 
mined through a set of coupled relations, each involving 
all the other t/-s, that have to be solved numerically. The 
corrections 5vj to the BCS pair-occupation amplitudes so 
produced further lower the ground state energy relative 
to projection after variation (but only by tenths of a per- 
cent). 

Extensive applications of such and related approaches 
in nuclear physics have led to the following conclusions: 
For reasonably small TV, as in nuclei, the explicit imple- 
mentation of projection techniques is tractable, though 
cumbersome. For very large TV they become intractable, 
but also unnecessary, since their corrections can be shown 
to vanish as TV -1 / 2 . However, even in nuclei the correc- 
tions to unprojected BCS theory are small (a few per- 
cent) in most cases, the only exception being very large 
couplings V > d. Thus, in most cases fixed-TV systems 
can perfectly adequately be described by BCS's grand- 
canonical wave function. Its TV-indefiniteness (and the 
associated breaking of gauge symmetry) then simply has 
the status of a clever calculational trick: it allows the use 
of a wave function so simple that the pair-occupation am- 
plitudes Vj can be found with a minimum of effort. The 
trick's justification is that the corrections Svj's produced 
by more careful approaches usually are small. (The de- 
vice of using symmetry-breaking wave- functions purely 
for the sake of calculational convenience is widespread 
in nuclear physics, and lucidly discussed in Ring and 
Schuck's bookEa in a chapter entitled "Restoration of 
Broken Symmetries".) 

The above conclusions imply that the following strat- 
egy should suffice for a qualitative description (more 
is not attempted here) of pairing correlations in iso- 
lated ultrasmall grains: although strictly speaking a 
fixed- TV technique would be appropriate, we shall adopt 
BCS's grand-canonical approach throughout, using Uj, Vj 
as grand-canonical approximations to Uj , Vj . Quantita- 
tively, this strategy is expected to become unreliable in 
the limit of large level spacing d/A > 1 (correspond- 
ing to "strong coupling" in nuclear applications). How- 
ever, the corrections due to a fixed-TV calculation (cur- 
rently under investigation applying projection!^ and ex- 
act diagonalizationEH methods), which should become 
significant in this regime, are not expected to be more 
severe than, for example, corrections arising from a 
non-equidistant level spectrum, which qualitatively are 
insignificantlij. 
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III. GENERALIZED VARIATIONAL BCS 
APPROACH 

Since in RBT's experiments T = 50mK <C d, A, we set 
T = 0. Our goal in this section is to calculate the discrete 
eigenenergies of an isolated, nm-scale metallic grain with 
pairing correlations, and understand their evolution in a 
magnetic field. To this end, we study the simplest con- 
ceivable pairing model within a generalized variational 
BCS approach. The results will be used in the next sec- 
tion as input into the calculation of the SET tunneling 
spectrum of such a grain (see Fig. |^ below). 



A. The Model 

The only symmetry expected to hold in realistic, 
irregularly-shaped ultrasmall grains at zero magnetic 
field is time-reversal symmetry. We therefore adopt a 
single-particle basis of pairs of time- reversed states |j±), 
whose discrete energies Sj are assumed to already incor- 
porate the effects of impurity scattering and the average 
of electron-electron interactions, etc. As simplest con- 
ceivable model describing a pairing interaction and a Zee- 
man coupling to a magnctic.finld, we adopt the following 
(reduced) BCS HamiltonianEEa: 



H = 



( £ 3-M+^)c] CT c iCT - \d^ i c) + c)_c r -c j >+ 



(6) 

Due to level repulsion the £j 's will be approximately uni- 
formly spaced. For simplicity, we take a completely uni- 
form spectrum with level spacing d, £j — jd+eo. Fluctu- 
ations in the level spacings have been studied with meth- 
ods of random matrix theoryli3, with qualitatively similar 
results. For a system with a total of N = 2m + p elec- 
trons, where the electron number parity p is for even TV 
and 1 for odd AT, we use the label j = for the first level 
whose occupation in the T — Fermi sea is not 2 but p. 

The pairing interaction is taken to include only states 
with \sj\ < uj c . Experimental evidence for the sufficiency 
of neglecting couplings between non-time-reversed pairs 
of states, i.e. of using only a reduced BCS-Hamiltonian, 
will be given in section IV D. For convenience we wrote 
the pair-coupling constant in Eq. (Q) as V — Xd, where 
A is a dimensionless parameter. The d — > "bulk gap" 
of the model thus is A = 2cj c e~ 1 / A . 

An applied magnetic field will completely penetrate 
an ultrasmall grain, since its radius (typically r ~ 5nm) 
is much smaller than the penetration length of 50 nm 
for bulk Al. The Zeeman term in Eq. (||), with ±h = 
±^iBgH, models the fact that the measured tunnel spec- 
tra of RBTEHa (shown in Fig. || in section IV B| ) evolve 
approximately linearly as a function of magnetic field, 
with g- factors between f .95 and 2 (determined from the 



differences between measured slopes of up- and down- 
moving lines). Deviations from g = 2 probably result 
from spin-orbit scattering, known to be small but non- 
zero in thin Al filmscd, but neglected below (where g = 2 
is used). Furthermore, orbital diamagnetism is also ncg=. 
ligible, just as for thin films in a parallel magnetic fieldEII 
but in marked contrast to bulk samples where it causes 
the Meissner effect: the grains are so small that even a 7T 
field produces a flux through the grain of only about 5% 
of a flux quantum O > which is too small to significantly 
affect the orbital motion of the electrons between sub- 
sequent reflections off the grain boundary. Some larger 
grains do show slight deviations from iJ-linearityo, which 
probably reflect .the onset of orbital magnetism (which 
gives correctionsES to the eigenenergies of the order of 
Hv Fr 3 (H I '0o) 2 ); however, these effects are much smaller 
than Zeeman energies in the grains of present interest, 
and will be neglected here. Thus, our model assumes that 
Pauli paramagnetism due to the Zeeman energy com- 
pletely dominates orbital diamagnetism, similarly to the 
case of thin films in parallel magnetic fieldsEj. 



B. The Variational Ansatz 

The Zeeman term favours states with non-zero total z- 
component of the total spin s = • s| (henceforth sim- 
ply called "spin"), so that increasing h will eventually 
lead to a series of ground state changes to states with 
successively larger spins. Therefore, we are interested 
in general in correlated states with non-zero spin, and in 
particular in their eigenenergies. We calculate these vari- 
ationally, using the following general Ansatz for a state 
|s, a) with a definite total, spin s (introduced by Soloviev 
for application in nucleiE2l): 



!">«> = IK« + II< 



(s,a) . (*,a) t t 
U) + Vf cj , ci_ 



|Vac). (7) 



The non-zero spin is achieved by placing 2s unpaired 
spin-up electrons in a set of 2s single particle states, 
say with labels j — a(l), a(2), ■ ■ ■ , a(2s) (see Fig. ||), 
while the remaining single-particle pairs of states have 
BCS-like amplitudes to be either filled (v^' 01 ^) or empty 

(u[ s ' a ^), with (u- s ' a ') 2 + (^ S:Q:) ) 2 = 1- The prime over 
products (and over sums below) indicates exclusion of 
the singly occupied states a(l), a(2), ■ ■ ■ , a(2s) (for which 
u (s,cx) ^ v (s,a) are n0 £ defined) , 

A short standard calculation reveals that the con- 
structed wave functions are orthogonal: (s,a\s',a') = 
Sss'Saa'- Therefore, the variational parameters v < ^ s ' a ' > 
and u^ s ' a ^ must be found independently for each (s, a) 
(hence the superscript). This is done by minimizing the 
variational "eigenenergies" 



£ s ,a(h,d) = (s,a\H\s,a) 
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| S! a)o = n c l (j -) + II c l c LlO). 



(14) 



i<0 



A^(«?' a) ) 4 



which we use to approximate the model's exact eigenen- 
ergies E 8t ct(h,d). Note that singly-occupied states are 
excluded from all primed sums involving uj's and Uj's. 
The last term, proportional to v 4 , is not extensive and 
hence neglected in the bulk case where only effects pro- 
portional to the system volume are of interest. Here we 
retain it, since in ultrasmall systems it is non-negligible 
(but not dominant either). 

Solving the energy-minimization conditions 

dS s>a /dv^ a) = (9) 
in standard BCS fashion yields 

(^ Q) ) 2 = (i-&/K| + < Q ] 1/2 )A (io) 

where the "pairing parameter" A s q< is determined by 
the generalized "gap equation" 



A \ .V^ ( s < a ) ( s - a ) 



A 2 



(11) 



(12) 



and £j = Ej — n — \d{ 



,(s,OL) 



Note that we retain the 



Xd{v 



(s,a)\2 



shift in usually neglected because it sim- 



ply renormalizes the bare energies, since for large d it 
somewhat increases the effective level spacing near Ep 
(and its neglect turns out to produce a significant upward 
shift in the £ s ,a(h, d)'s, which one is trying to minimize). 
The chemical potential /i is fixed by requiring that 



2m 



P = (a,a\N\s,a) = 2s + 2^(^ a) ) 2 



(13) 



Generally Eqs. ©, @) and @ have to be solved si- 
multaneously numerically. In the limit d/A — > (inves- 
tigated analytically in Appendix A 1 ) , Eq. ( |l2| ) reduces 
to the standard bulk T = gap equation. 

In contrast to conventional BCS theory, the pairing 
parameter A S! ct can in general not be interpreted as an 
energy gap and is not an observable. It should be viewed 
simply as a mathematical auxiliary quantity which was 
introduced to conveniently solve Eq. (^|). However, by pa- 
rameterizing the variational quantities v^ s,OL ^ and u^ B,a \ 
A s jy, does serve as a measure of the pairing correlations 
present in a state |s, a), since for vanishing A s a the lat- 
ter reduces to an uncorrelated paramagnetic state with 
spin s, namely 



We shall denote the energy of this uncorrelated state by 
f s ° Q = o(s,ol\H\s,ol)o, and define the "correlation en- 
ergy" of \s, a) as the energy difference £fa = £ s ,a — 

fO 



C. Qualitative Discussion 

Before launching into numerical results, let us antici- 
pate by qualitative arguments what is to be expected: 

Firstly, the gap equation for A Si a{d) is h- independent. 
The reason is that only those j-levels contribute in the 
gap equation that involve correlated pairs of states, each 
of which have spin and hence no Zeeman energy. Con- 
sequently, the — 2s/i-dependence of £ s ,a in Eq. || is simply 
that of the 2s unpaired electrons. 

Secondly, the discreteness of the sum in the gap equa- 
tion ( |l2|) will cause A Sy ct to decrease with increasing d. 
To see this, inspect Fig. [jj in which the height of each ver- 
tical line represents the value of UjVj for a time-reversed 
pair |j±). Figs. |I](a) to (c) illustrate that an increase in 
level spacing implies a decrease in the number of pairs 
with significant pair- mixing, i.e. those within A of £f 
which have non-zero UjVj. This number can roughly 
speaking be called the "number of Cooper pairs" of the 
system. Since for d ^> A no pairs lie in the correlated 
regime \sj — £p\ < A where pair-mixing occurs, A St a 
will be zero in this limit, so that in general A S: a(d) will 
be a decreasing function of d, dropping to zero at about 
d ~ A. Physically speaking, this happens since with 
increasing d the increasing kinetic energy cost of pair- 
mixing (which shifts pair-occupation amplitude from be- 
low to above £f) causes the correlations to weaken, be- 
coming negligible for large enough d. 

Thirdly, the (s, o:)-dependent restriction on the primed 
sum in the gap equation implies that A s- a{d) at fixed d 
will decrease with increasing s: larger s means more un- 
paired electrons, more terms missing from the primed 
sum, less correlated pairs and hence smaller A St a- Tha 
physics behind this has been called the blocking effecvEl 
in nuclear physics: Singly-occupied states cannot take 
part in the pair-scattering caused by the BCS-like inter- 
action (^) and hence decrease the pha se space for pair 
scattering, as explained in section II A. (Their absence 
in the primed sum simply reflects this fact.) The block- 
ing effect becomes stronger with increasing d, since then 
the relative weight of each term missing in the primed 
sum increases. It also is stronger the closer the blocked 
state lies to £f, since the excluded Uj' a ^Vj' a ^ contribu- 
tion to theprimed sum is largest near ef, as is evident 
from Fig. [j]. On the other hand, an unpaired electron 
will have almost no blocking effect if \£j — Ep\ ^ A, since 



u (s,cx) v (s,a) ygjjjgjjgg there anyway. 
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Finally, note that the (s, a)-dependence of A s>a for 
d ~ A illustrates why in this regime a conventional mean- 
field treatment is no longer sufficient: the system cannot 
be characterized by a single pairing parameter, since the 
amount of pairing correlations vary from state to state, 
each of which is characterized by its own pairing param- 
eter. 



D. General Numerical Solution 

It is possible to solve the modified gap equation an- 
alytically in two limits, d <C A and d 3> A s (see Ap- 
pendix [A]) , but generally the gap equation and (|l3|) have 
to be solved numerically. In doing so, some assump- 
tions are necessary about parameter values (though us- 
ing slightly different values would not change the results 
qualitatively). We measure all energies in units of the 
bulk gap A = 2oj c e~ 1 / x of the model. However, its ex- 
perimental value differs from that of a truly bulk-,system, 
since it is known from work with Al thin filmsQiHa that 
the effective dimensionless pairing-interaction strength A 
is larger in Al samples of reduced dimensionality than 
in truly bulk, three-dimensional systems. (Though true 
for Al, this is not a universal property of small samples, 
though — for Nb, A is larger in the bulk than in thin 
filmsE-3.) Since thin films in a parallel magnetic field are 
analogous in many ways to ultrasmall grains (in particu- 
lar regarding the dominance of Pauli paramagnetism over 
orbital diamagnetism) , we shall assume that the effective 
coupling constant A is the same in both. Adopting, there- 
fore, the value A = 0.38meV found for thin Al films in 
Ref. |H, and taking the cut-off to be the Debye frequency 
uj c = 33meV of Al, we use A = [\n{2u> c / A)]- 1 = 0.224 for 
the dimensionless pairing-interaction strength. Further- 
more, we smeared the cutoff of the BCS interaction over 
two single-electron levels, to ensure that discontinuities 
do not occur in d-dependent quantities each time the en- 
ergy \ej = dj + Eq\ of some large-|j| level moves beyond 
the cut-off uj c as d is increased. 

Solving Eqs. ©, (|TJ|) and dfj) is a straightforward 
numerical exercise which we performed, for the sake 
of "numerical consistency" , without further approxima- 
tions. (Since some minor approximations were made in 
Ref. |f5[ e.g. dropping the Acfo| term in £j, and slightly 
different parameter-values were used, the numerical re- 
sults there sometimes differ slightly from the present 
ones; see, e.g. Fig. ||.) It should be understood, though, 
that only qualitative significance can be attached to our 
numerical results, since our model is very crude: it ne- 
glects, for instance, fluctuations in level spacing and in 
pair-coupling constants, and we do not carry out a fixed- 
N projection, all of which presumable would somewhat 
influence the results quantitatively. 



J]K + <cJ + c]_)|Vac). 



(15) 



1. Spins Ground States 

In a given spin-s sector of Hilbert space (with p = 
2s mod 2), let s) be the variational state with the lowest 
energy, i.e. the "variational spin-s ground state" . It is 
obtained by placing the 2s unpaired electrons as close 
as possible to £f [Fig. ||(a)], because this minimizes the 
kinetic energy cost of having more spin ups than downs: 

s-l+p/2 

w = n «3h 

j=-s+p/2 % 

(The particular choice of a. in the general Ansatz (Q) to 
which |s) corresponds is a(n) = n— [s] — I for n = 1 . . . 2s, 
where [s] is the largest integer < s.) The numerical 
results for the corresponding pairing parameters A s (d), 
shown in Fig. |3](a) for some several small s, confirm the 
properties anticipated in the previous subsection's qual- 
itative discussion: 

Firstly, each A s decreases with d, vanishing at a crit- 
ical level spacing d CjS beyond which no pair-mixing cor- 
relations exist in this level of approximation. In Ap- 
pendix A 2 it is shown that near d c , s , A s (d) has the 
standard mean-field form yl — d/d CtS ; this was to be 
expected, since the variational approach to finding |s) is 
equivalent to doing standard mean-field theory within the 
spin-s sector of Hilbert space. (Note that one should not 
attach too much significance to the precise numerical val- 
ues of the d c , s reported in Fig. ||, since they depend sen- 
sitively on model assumptions: for example, the values 
for d c s) and <i c ,i/2 differ somewhat from those reported 



in Refs. |f2|,[l5|, due to their use of a slightly different A 
and minor numerical approximations not used here, as. 
mentioned above. Moreover, Smith and Ambegaokaro 
showed that the precise distribution of levels used influ- 
ences d c , s significantly.) 

Secondly, A s decreases rapidly with increasing s at 
fixed d (and d c , s < c? c , S ' if s > s'), illustrating the block- 
ing effect. This result, which is expected to be indepen- 
dent of model details, is a generalization of the parity 
effect discussed by von Delft et alS3. (They studied only 
ground state pairing correlations and found that these 
are weaker in odd (s = 1/2) grains than in even (s = 0) 
grains, A odd = A 1/2 < A even = A .) The blocking ef- 
fect is most dramatic in the regime d/A € [0.77,2.36] 
in which Ao ^ but A s ^oi-t= 0- This is a regime of 
"minimal superconductivity" E3, in the sense that all pair- 
ing correlations that still exist in the even ground state 
(since Ao ^ 0) are completely destroyed by the addition 
of a single electron or the flipping of a single spin (since 
A s ^o = 0). 

Fig. ||(b) shows the eigenenergies £ s (solid lines) of |s) 
and the energies £° (dotted lines) of the corresponding 
uncorrelated paramagnetic states 



-l+p/2 



i s >°= n 4+ n c i c Livac> 
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The solid and dashed spin-s lines meet at the critical level 
spacing c? CjS , above which no pairing correlations survive. 



2. Spin-s Excited States 

Among all possible excited states with definite s, we 
consider here only those created from |s) by exiting one 
electron from the topmost occupied level s — 1 + p/2 of 
s) to some higher level j + s — 1 + p/2: 



\*,J) 



s-2+p/2 

nc f 
] + 

] = -s+p/2 

I 

xn«+^M+ c l-)|Vac) 



(17) 
(18) 



(This reduces to \s) if j = 0; the particular choice of a. in 
Ansatz (0) to which \s,j) corresponds is a (n) = n—[s] — l 
for n = 1 ... 2s - 1 and a(2s) = [s] - 1 + j.) 

Interestingly, one finds that the larger j, the longer the 
pairing correlations survive with increasing d. This is il- 
lustrated by the simple example s = 1/2: Fig. |^(a) shows 
that the critical spacings d c i/ 2 .j (at which the pairing pa- 
rameters Ai/2, j(d) vanish) increase with j, approaching 
the value d c .o of the spin-0 case as j — > oo. This re- 
sult is reflected in the excitation energies of Fig. ||(b): 
the excited states of a the spin-1/2 sector have non-zero 
correlation energies (difference between solid and dashed 
lines) at d- values for which the spin-1/2 ground state cor- 
relation energy of Fig. ||(b) is already zero. The intuitive 
reason why more-highly-excited states have more pair- 
ing correlations than the corresponding spin-1/2 ground 
state 1 1/2) is of course quite simple: The larger j, i.e. the 
further the unpaired electron sits from the Fermi sur- 
face where pairing correlations are strongest, the less it 
disrupts pair- mixing (since UjVj becomes very small for 
large j, see Fig. |l|). In fact, for very large j, the state 
will have the just about same amount of pairing 
correlations as the even ground state |0) (A^j — Ao), 
since the unpaired electron sits so far from ef that the 
pairing correlations are effectively identical to those of 

Similar effects are seen for excited states in other spin 
sectors s ^ |. The higher the excitation, the larger the 
pairing parameter A Sy ot- Nevertheless the energy of the 
excited states is always higher than that of the corre- 
sponding spin-s ground state, since the kinetic-energy 
cost of having an unpaired electron far from sf can be 
shown to always outweigh the interaction-energy gain due 
to having less blocking and hence a larger A s a . 



E. Magnetic Field Behaviour 

In a magnetic field, the Zeeman energy favors states 
with non-zero spin. However, since such states have 



smaller correlation energy due to the blocking effect a 
competition arises between Zeeman energy and correla- 
tion energy. The manifestations of the blocking effect 
can thus be probed by turning on a magnetic field; if it 
becomes large enough to enforce a large spin, excessive 
blocking will destroy all pairing correlations. 

The situation is analogous to ultra-thin films in a 
parallel magnetic field,EiJ where orbital diamagnetism is 
negligible for geometrical reasons and superconductiv- 
ity is destroyed at sufficiently large h by Pauli para- 
magnetism. This occurs via a first order transition 
to a paramagnetic state, as predicted by Clogston and 
Chandrasekhar (CC)EJt!l by the following argument (for 
bulk systems): A pure Pauli paramagnet has ground 
state energy —h 2 Af(ep) and spin s = hJV(ep) (since it 
chooses its spin such that the sum of the kinetic and Zee- 
man energies at spin s, s 2 N{ef) — 2hs, is minimized). 
When this energy drops below the bulk correlation en- 
ergy — ^A 2 Af(ep) of the superconducting ground state, 
which happens at the critical field hcc = A/ y/2, a transi- 
tion will occur from the superconducting to the paramag- 
netic ground state. The transition is first-order, since the 
change in spin, from to sec = h cc M{e F ) = A/(dV2), 
is macroscopically large (AT(sf) = l/d ~ Vol). In 
tunneling experiments into ultra-thin (5nm) Al films 
(A = 0.38mc)Land H cc = 4.7T) this transition has 
been observedEj as a jump in the tunneling threshold 
(from A — hcc to zero) at hcc- 

In isolated ultrasmall grains, the above picture of the 
transition needs to be rethought in two respects due to 
the discreteness of the electronic spectrum: Firstly, the 
spin must be treated as a discrete (instead of continu- 
ous) variable, whose changes with increasing h can only 
take on (parity-conserving) integer values. Secondly, one 
needs to consider more carefully the possibility of h- 
induced transitions to non-zero spin states that are still 
pair- correlated (instead of being purely paramagnetic), 
such as the variational states \s, a) discussed above. (In 
the bulk case, it is obvious that such states play no role: 
the lowest pair-correlated state with non-zero spin ob- 
tainable from the ground state by spin flips is a two-quasi- 
particle state, costing energy 2 A— 2ft,; when h is increased 
from 0, the paramagnetic transition at hcc = A/v2 thus 
occurs before a transition to this state, which would re- 
quire h = A, can occur.) 

Within our variational approach, the effect of increas- 
ing h from can be analyzed as follows: At given d and h, 
the grain's ground state is the lowest-energy state among 
all possible spin-s ground states |s) having the correct 
parity 2s mod 2 = p. Since £ s (h,d) = £ s (0,d) — 2hs, 
level crossings occur with increasing h, with £ s i dropping 
below £ s at the level crossing field 



.(d) 



£ s ,{0,d)-£ s (0,d) 
2(s'-s) 



(19) 



Therefore, as h is slowly turned on from zero with ini- 
tial ground state \s = p/2), a cascade of successive 
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ground-state changes (GSC) to new ground states |si), 



*2 



.will occur at the fields h SOtS1 , h si> 



. We de- 



note this cascade by (sq, si); (si, S2); . . ., and for each of 
its ground state changes the corresponding level-crossing 
fields h StS >(d) is shown in Fig. |[ Generalizing CC's criti- 
cal field to non-zero d, we denote the (parity-dependent) 
field at which the first transition (so,si) occurs by 
hcc(d,p) = h SOtS1 (d), which simply is the lower envelope 
of the level-crossing fields h So , sl in Fig- ||- In the limit 
d — ► wc find numerically that it correctly reduces to 
the Clogston-Chandrasekhar value hcc(0,p) — A/y/2. 

In general, the order in which the GSCs occur with in- 
creasing h depends sensitively on d and an infinite num- 
ber of distinct regimes (cascades) I, II, III, . . . can be 
distinguished: Starting at large d we find the typical 
normal behaviour (0, 1); (1, 2); (2, 3); . . . for even grains 
and (|, |); (|, |); . . . for odd grains, with ho t i < (or 
>) I in regimes I (or II). In regimes III and IV of 
somewhat smaller d, the order of GSCs is (0, 2); (2, 3); . . . 
and (|, |); (|, |); . . ., etc, i.e. the spin si attained af- 
ter the first GSC (soi s i) has increased to 2 in the even 
case. This illustrates a general trend: the spin si(d) af- 
ter the first transition increases with decreasing d and 
becomes macroscopically large in the d — > limit, where 
si = hcc/d = A/ (d\/2), as explained in recounting CC's 
argument above. 

Furthermore, it turns out that A Sl (d) = for all d, 
implying that after the first GSC the new ground state 
\si) is always (not only in CC's bulk limit) an uncorre- 
lated, purely paramagnetic state. In this regard, CC's 
picture of the transition remains valid throughout as d 
is increased: at hcc(d,p), a transition occurs from the 
superconducting ground state to a paramagnetic, uncor- 
related state |si)o, the transition being first-order in the 
sense that A Sl (d) = 0; however, the first-order transition 
is "softened" with increasing d, in the sense that the size 
of the spin change, s% — sq, decreases from being macro- 
scopically large in the bulk to being equal 1 at d 3> A 
(regimes I and II). 



F. Deficiencies of the Variational Ansatz 

Though the variational method we used to calculate 
the systems "eigenenergies" is expected to yield qualita- 
tively correct results, it does have some deficiencies: 

Firstly, a variational approach by construction only 
gives an upper bound on the exact eigenenergies -E^a- 
The variational energies £ s ,a could be lowered further by 
choosing better trial wave functions that sample larger 
parts of a given spin-s Hilbert space, i.e. by including 
"fluctuations" about the chosen states. 

Secondly, the abrupt vanishing of the pairing param- 
eters A s ,a(d) ~ \J\ — d/d CtS at a critical level spacing 
d s ,a. [see Appendix A 2 and Fig. ||(a))] is unphysical: in 
a finite system, any non-zero pair-interaction constant 
will always induce a non-zero amount of pairing corre- 



lations, i.e. the canonical A' s {d) of Eq. (Q) will always 
be non-zero, though it could become arbitrarily small for 
sufficiently large d. (This statement is analogous to stat- 
ing that "in a finite system no abrupt phase transition 
between a zero and non-zero order parameter occurs."). 
The abrupt, mean-field-like vanishing of A St a(d) is of 
course an artefact, that occurs since the grand-canonical 
variational Ansatz is equivalent (at least for the spin-s 
ground states |s)) to doing mean- field theory in a fixed- s 
Hilbert space. 

Thirdly, the variational states of course are not N- 
eigenstates (though they do have definite parity), and 
Eq. ( |lg| ) only fixes the mean electron number. Our 
reasons for nevertheless adopting t hem to describe an 
isolated grain were given in section [I C : a large body 
of experience in nuclear physics showed that fixed-TV 
projections generally produce only minor corrections to 
the grand-canonical BCS results. Nonetheless, note 
that we expect a fixed-A projection (currently under 



investigationEj) to somewhat ameliorate the first two of 
the above-mentioned deficiencies of the variational ap- 
proach: projection after variation of \s) to fixed N will 
lower the energy £ s a bit, and presumably projection be- 
fore variation will in addition result in a canonical pair- 
ing parameter A s (g0 that decays smoothly with increas- 
ing d from finite to arbitrarily small but non-zero values. 
Note, though, that this is not expected to change the 
eigenenergies very much, since the correlation energies 
rapidly approach zero anyway when the correlations be- 
come weak. In other words, we expect the variational 
scheme for calculating eigenenergies to break down only 
when A s becomes so small that it has no experimen- 
tal relevance any more (to check_this in detail, strictly 
canonical calculations are neededEjO) . 



IV. OBSERVABLE QUANTITIES 

In this section, we consider the grain coupled to leads 
as in RBT's SET experiments. After explaining what 
kind of information can and can not be extracted from 
their data, we turn to the calculation of observable quan- 
tities: (a) We calculate theoretical tunneling spectra and 
compare these to RBT's measurements; (b) address the 
question of the observability of various parity effects, 
proposing to search for one involving the pair-breaking 
energy; and (c) explain how RBT's data give direct evi- 
dence for the dominance of time-reversed pairing, at least 
for small fields. 



A. Experimental Details 

In RBT's experimentsiflj§, an ultrasmall grain was 
used as central island in a SET: it was connected via 
tunnel barriers to external leads and capacitively cou- 
pled to a gate, and its electronic spectrum determined 
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by measuring the tunnel current through the grain as a 
function of transport voltage (V), gate voltage (V g ) and 
magnetic field (H = Ii/hb, with hb = 0.0571meV/T) at 
a fixed temperature of 50mK. 

The particular grain (Ref. |J Figs. 1(b) ,2,3) with which 
we shall compare our theory had the following param- 
eters: Its radius was estimated as r ~ 4.5nm by as- 
suming the grain to be hemispherical, implying a vol- 
ume ~ (5.7nm) 3 and a total number of conduction elec- 
trons N of about 3 x 10 4 . The crude order-of-magnitude 
free-electron estimate d = 2-k 2 T% 2 / (mfcpVol) for the mean 
level spacing near £p yields d ~ 0.45meV. The SET had 
lead-to-grain capacitances C\ = 3.5aF, C2 = 9.4aF, gate- 
to-grain capacitance C g — 0.09aF and charging energy 
Ec = e 2 /2Ctotai = 46meV. The tunnel current is on the 
order of 10 _10 A, implying an average time of 2 • 10 _9 sec 
between subsequent tunneling processes. 

Since the charging energy Ec was very much larger 
than all other energy scales, such as the bulk gap 
(A ~ 0.38meV), typical values of the transport voltage 
(V~lmV) and the temperature, fluctuations in electron 
number on the grain are strongly suppressed, so that co- 
herent superpositions between states with different N 
need not be considered. The energy-balance condition 
that determines through which eigenstates of the grain 
electrons can tunnel for given values of transport- and 
gate voltage thus involve differences between the eigenen- 
ergies of a grain with fixed particle number N or N ± 1 , 

AE if = {Ef + E^)- (£f ±J + E^ 1 ) , (20) 

corresponding to the energy cost needed for some rate- 
limiting electron tunneling process |z) jv±i - * |/)jv off or 
onto the grain. Here |/)jv denotes a discrete eigenstate 
of the ./V-electron grain with eigenenergy Ey 7 + Eq . Fol- 
lowing the "orthodox model" of SET charging, we take 
Eq , the grain's electrostatic energy (relative to a neutral 
grain with iV electrons) as Eq = E C (N — N — Q g /e) 2 , 
where Q g — C g V g + const is the gate charge, and assume 
the Coulomb-interaction to be screened sufficiently well 
that its sole effect is to shift all fixed-iV eigenstates by 
the same constant amount Eq. (The latter assumption 
is somewhat precarious: it becomes worse with decreas- 
ing grain size, and was shown to break down in grains 
half the present size. El) 

RBT were able to extract the energy differences AE{ f 
from their data: the differential conductance dl/dV as 
function of V at fixed V g has a peak whenever eV times 
a known capacitance ratio is equal to one of the A£^'s, 
at which point another channel for carrying tunneling 
current through the grain opens up (the inclusion of the 
capacitance ratio takes into account that the voltage drop 
across each of the two tunnel junctions can be different 
if their capacitances are not identical]) . Plotting the po- 
sition of each conductance peak as function of h gives 
the so-called experimental tunneling spectrum shown in 
Fig. ^, in which each line reflects the ^-dependence of 
one of the energy differences AE(h). 



It is important to note that the experimental threshold 
energy at ft, = for the lowest-energy tunneling process 
(y-intercept of the lowest line, the so-called "tunneling 
threshold" ) yields no significant information, since it de- 
pends on the grain's change in overall charging energy 
due to tunneling, 

6E C =Eg- E^ 1 = E c [Q g /e ~(N-N ± ±)], (21) 

which depends (via Q g ) in an imprecisely-known way on 
the adjustable gate voltage V g . This V^-dependence can 
usually (e.g. in SETs with much smaller charging ener- 
gies than here) be quantified precisely by studying the 
Coulomb oscillations that occur as function of V g at fixed 
V. Unfortunately, in the present case a complication 
arisesc-a due to the smallness of the gate capacitance: to 
sweep Q g through one period of 2e, the gate voltage V g 
must be swept through a range so large (2e/C g ~ 3.5V) 
that during the sweep small "rigid" shifts of the en- 
tire tunneling spectrum occur at random values of V g . 
They presumably are due to single-electron changes in 
the charge contained in other metal grains in the neigh- 
borhood of the grain of interest; these changes produce 
sudden shifts in the electrostatic potential of the grain, 
and thus spoil the exact 2e-periodicity that would other- 
wise have been expected for the spectra. 

In contrast to the threshold energy, however, the sep- 
arations between lines, 

AEt f , - AE if = Ef, - Ef , (22) 

are independent of gate voltage and hence known abso- 
lutely; they simply correspond to the differences between 
eigenenergies of a fixed-N grain, i.e. give its fixed-N ex- 
citation spectrum, and these are the quantities that we 
shall focus on calculating below. 

The most notable feature of RBT 's measured tunneling 
spectra is the presence (absence) of a clear spectroscopic 
gap 2fl e > d between the lowest two lines of the odd- 
to-even (even-to-odd) measured spectra in Fig. ||(a,b). 
This reveals the presence of pairing correlations: in even 
grains, all excited states involve at least two BCS quasi- 
particles and hence lie significantly above the ground 
state, whereas odd grains always have at least one quasi- 
particle and excitations need not overcome an extra gap. 

Since the AE^f's in Eq. ( pc| ) are constructed from 
fixed-N and fixed-N ± 1 eigenenergies, we shall approx- 
imate these using the variational energies £ s a discussed 
in previous sections for a completely isolated grain. (We 
thereby make the implicit assumption that the grain's 
coupling to the leads is sufficiently weak that this does 
not affect its eigenenergies, i.e. that the leads act as 
"ideal" probes of the grain.) The £ s a will be used as 
starting point to discuss various observable quantities; in 
particular, we shall make contact with RBT's experimen- 
tal results by constructing the theoretical tunnel spec- 
trum (as function of h and d) predicted by our model. 
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B. The Tunneling Spectrum in an Magnetic Field 

The kind of tunneling spectrum that results depends 
in a distinct way on the specific choice of level spacing 
d and final-state parity p (i.e. the parity of the grain 
after the rate- limiting tunneling process has occured). 
To calculate the spectrum for given d and p, we proceed 
as follows below: we first analyze at each magnetic field 
h which tunneling processes |i)jv±i ~~ * \/)n are possible, 
then calculate the corresponding energy costs AEif(h) 
of Eq. (po|) and plot AEij(h) — AE mm (0) as functions 
of h for various combinations of i, /, each of which gives 
a line in the spectrum. We subtract AE m i n (0), the h = 

threshold energy cost for the lowest-lying transition, 
since in experiment it depends on V g and hence yields no 
significant information, as explained above. 

Fig. [7] shows four typical examples of such theoretical 
tunneling spectra, with some lines labeled by the corre- 
sponding \i) — ► I/) transition. 

When taking the data for Fig. |], RBT took care to 
adjust the gate voltage V g such as to minimize non- 
equilibrium effects, which we shall therefore neglect. For 
given h, we thus consider only those tunneling processes 
for which the initial state \i) corresponds to the grain's 
ground state |sj) at that h (and d,p), whose spin Sj can 
be inferred from Fig. 0. Since the grain's large charg- 
ing energy ensures that only one electron can tunnel at 
a time, the set {|/)} of possible final states satisfies the 
"spin selection rule" \sf — s,| = i and includes, besides 
the spin-s/ ground state |s/), also excited spin-s/ states. 

Whenever h passes through one of the level-crossing 
fields h Si _ s ., of (119), the grain experiences a ground state 
change (sj,Sj'). After this GSC, |sj>) is the new initial 
state for a new set of allowed tunneling transitions \si>) — » 
{I s /')} (satisfying I s/' — I = 1/2). Since this new set in 
general differs from the previous set of transitions |sj) — > 
{|/)} allowed before the GSC, at h SitS ., one set of lines in 
the tunneling spectrum ends and another begins. A line 
from the former connects continuously to one from the 
latter only if its final state |/) can be reached from both 

1 Si) and I s^) [i.e. if Sf — Si — — (sf — Sj')]; in this case, 
the two lines |sj) — > |/) and |sj') — * |/) join at h sus '. 
via a kink, since AEif(h) and AEif(h) have slopes of 
opposite sign. However, for most lines this is not the 
case (since usually |s/ — Sj/| ^ 1/2), so that at h Si>s > the 
line I Si) — > \f) simply ends while new lines |si') \f) 
begin. This results in discontinuities (or "jumps" ) in the 
spectrum at h SitS { of size {AEvfi — AEif)(h Si s <), unless 
by chance some other final state \f) happens to exist for 
which this difference equals zero. 

Since the order in which the GSCs (s^s-j') occur as 
functions of increasing h depend on d and p, as indicated 
by the distinct regimes I, II, III, ... in Fig. ||, one finds 
a distinct kind of tunneling spectrum for each regime, 
differing from the others in the positions of its jumps 
and kinks. In regime I, where the order of occurrence of 
GSCs with increasing h is (0, 1); (|, §); (1,2); (f , §); . . ., 



there are no discontinuities in the evolution of the lowest 
line [see Fig. 0(a)]. For example, for the e-*o spectrum, 
the lowest |0) — > |l/2) line changes continuously to 1 1) — > 
|l/2) at /io,ij si nce I s / — s il = 1/2- However, in all other 
regimes the first change in ground state spin (at ho tSl 
from to si) is > 1, implying a jump (though possibly 
small) in all e — > o lines, as illustrated by Fig. 0(b). 

The jump's magnitude for the tunneling thresholds, 
i.e. the lowest e — > o and o — > e lines, is shown as func- 
tion of d in Fig. M. It starts at d = from the CC 
value A(l — l/\/2) measured for thin Al films,u and 
with increasing d decreases to (non-monotonically, due 
to the discrete spectrum). This decrease of the size of 
the jump in the tunne ling threshold reflects the fact, dis- 
cussed in section [II E , that the change in spin at the first 
ground state change (so,si) decreases with increasing d 
(as Si — so ~ hec/d), and signals the softening of the 
first-order superconducting-to-paramagnetic transition. 

The fact that the measured tunneling thresholds in 
Fig. |^ show no jumps at all (which might at first seem sur- 
prising when contrasted to the threshold jumps seen at 
hec m thin films io-a parallel fielder) , can therefore nat- 
urally be explainedt-3 by assuming the grain to lie in the 
"minimal superconductivity" regime I of Fig. (where 
the jump size predicted in Fig. @ is zero). Indeed, the 
overall evolution (i.e. order and position of kinks, etc.) of 
the lowest lines of Fig. || qualitatively agrees with those 
of a regime I tunneling spectrum, Fig. [7|(a). This allows 
us to deduce the following values for the level-crossing 
fields H s . s ' (indicated by vertical dashed lines in Figs, g 
and @): Hq A = 4T, H 1/2 , 3/2 = 4.25T, H ia = 5.25T and 
^3/2.5/2 = 6.5T. As corresponding uncertainties we take 
AH S ' US >. = 0.13T, which is half the H resolution of 0.25T 
used in experiment. 

By combining the above H SuS > values with Fig. ||, 
some of the grain's less-well-known parameters can be 
determined somewhat more precisely: Firstly, the grain's 
"bulk Hcc" neld can be estimated by noting from Fig. || 
that h ,i/h C c - 0.95, so that H cc = ff ,i/0.95 ~ 4.2T. 
This is in rough agreement with the value Hcc — 4.7T 
found experimentallycj in thin films in a parallel field, 
confirming our expectation that these correspond to the 
"bulk limit" of ultrasmall grains as far as paramag- 
netism is concerned. (Reca ll that our numerical choice of 
A = 0.224 in Section [II D was based on this correspon- 
dence.) Secondly, the grain's corresponding bulk gap is 
A = V^I-LbHcc - 0.34meV. Thirdly, to estimate the 
level spacing d, note that since H1/2, 3/2/^0,1 — 1-06, 
this grain lies just to the right of the boundary between 
regions II and I in Fig. || where d/A ~ 0.63, i.e. d ~ 
0.21meV. ( The c rude volume-based value d ~ 0.45meV 
of Section [V A| thus seems to have been an overesti- 
mate.) It would be useful if the above determination 
of d could be checked via an independent accurate ex- 
perimental determination of d directly from the spacing 
of lines in the tunnel spectrum; unfortunately, this is not 
possible: the measured levels are shifted together by in- 
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teractions, implying that their spacing does not reflect 
the mean independent-electron level spacing d. 

The higher lines plotted in Fig. correspond to states 
where the electron tunnels into an excited spin-s f state. 
For simplicity we considered only excited states \sf, j) in- 
volving a single electron- hole excitation rel ative to \sf), 



1. Exp — Eo is not measurable at present 



or as 



such as the example discussed in section III D 2 
sketched in Fig. ||(b), though in general others are ex- 
pected to occur too. The jumps in these lines (e.g. in 
Fig. |^(a) at ft.1,2) occur whenever the two final excited 
states \sf,jf) and \sf>,jf>) before and after the GSC at 
h s . s > have different correlation energies. (Recall that 
the correlation energy of an excited state |s/,ct/) can 
be non-zero even if that of the corresponding ground 
state \sf) is zero, since the former's unpaired electrons 
are furt her aw ay from ep, so that A Sfi a. f > A s ,, see 
Section III D 2| .) Experimentally, these jumps have not 
been observed. This may be because up-moving, reso- 
nances lose amplitude and are difficult to followtJ with 
increasing h, or because the widths of the-.excited reso- 
nances (~ 0.13 A) limit energy resolution.lHJ 

For somewhat larger grains, the present theory pre- 
dicts jumps even in the lowest line. It remains to be in- 
vestigated, though, whether orbital effects, which rapidly 
increase with the grain size, would not smooth out such 
jumps. 

Finally, note that more than qualitative agreement be- 
tween theory and experiment can not be expected: in 
addition to the caveats mentioned in the second para- 
graph of Section III D[ , we furthermore neglected non- 
equilibrium effects in the tunneling process and assumed 
equal tunneling matrix elements for all processes. In re- 
ality, though, random variations of tunneling matrix el- 
ements could suppress some tunneling processes which 
would otherwise be expected theoretically. 



One might expect that the odd-even ground state en- 
ergy difference Eq B = {E 1 / 2 — Eq) should reveal traces of 
the parity effect. Unfortunately, at present this quantity 
is not directly measurable, for the following reasons: 

If the transport voltage V is varied at fixed gate- 
voltage V g , the energy cost of changing the grain's elec- 
tron number by one (the h — threshold tunneling en- 
ergy) depends [see Eq. (|20|)1 not only on Eq C but also on 
the change SEc in the grain's charging ene rgy d ue to tun- 
neling. However, as explained in Section IV A, SEc de- 



pends (in an imprecisely-known way) on the actual value 
of V g . Therefore only the grain's fixed- N excitation spec- 
trum (distance between lines of tunneling spectrum) can 

be measured accurately in this way, but not Eq £ . 

If the gate voltage V g is varied at a fixed transport volt- 
age in the linear response regime V ~ 0, i.e. Coulomb os- 
cillations are studied, one expects to find a 2e-periodicity 
in the so-called gate charge Q g = C g V g + const, with 

Eq £ determining the amount of deviation from the 
e-periodicity (see Appendix [b| for details). This— lias 
been demonstrated convincingly in ^m-scale devicestElEl 
However, for nm-scale devices it is at present not possible 
to study (as suggested in Ref. |lj) e— or 2e-dependent fea- 
tures with sufficient accuracy to carry through this pro- 
cedure: the charging energy is so large that the pairing- 
induced deviation from e-periodicity is a very small effect 
(a fractional change of order Eq / Ec < 0.01), which can 
easily be obscured by V^-dependent shifts in backg round 
charge near the transistor, as explained in Section IV A. 



Parity Effect in Pair-Breaking Energies 



C. Parity Effects 

As mentioned in the introduction, several authorstHrt^l 
have discussed the occurrence of a parity effect in ul- 
trasmall grains: "superconductivity" (more precisely, 
ground state pairing correlations) disappears sooner with 
decreasing grain size in an odd than an even grain 
(A^ < Ao, and d c i/ 2 < g? c ,o)- This is a consequence of 
the blocking effect, which is always stronger in the pres- 
ence of an odd, unpaired electron than without it. This 
section is devoted to discussing to what extent this and 
related parity effects are measurable. 

Since pairing parameters such as A]y 2 , Ao are not ob- 
servable quantities, measurable consequences of parity 
effects must be sought in differences between eigenener- 
gies, which are measurable. 



Since the quantities that RBT can measure accurately 
are fixed-iV excitation spectra, let us investigate what 
parity-effects can be extracted from these. Since any par- 
ity effect is a consequence of the blocking effect, we begin 
by discussing the latter's most obvious manifestation: it 
is simply the fact that breaking a pair costs correlation 
energy, since the resulting two unpaired electrons disrupt 
pairing correlations. This, of course, is already incorpo- 
rated in mean-field BCS theory via the excitation energy 
of at least 2A involved in creating two quasi-particles. 
It directly manifests itself in the qualitative difference 
between RBT's even and an odd excitation spectra (ex- 
plained in section IV A), namely that the former shows 
a large spectral gap 2£l e > d between its lowest two lines 
that is absent for the latter (Fig. ^|). ■— ■ 

The parity effect discussed by von Delft et alS3 and 
Smith and AmbegaokarH3 referred to a more subtle con- 
sequence of the blocking effect that goes beyond conven- 
tional BCS theory, namely that the pairing parameters 
A s have a significant s-dependence once d/A becomes 
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sufficiently large. Although these authors only consid- 
ered the ground state parity effect Ai/ 2 < A , the same 
blocking physics will of course also be manifest in gen- 
eralizations to s > h. In fact, the problems with mea- 
suring the odd-even ground state energy difference Eq C 
discussed above leave us no choice but to turn to s > i 
cases when looking for a measurable parity effect. Specifi- 
cally, we shall now show that a parity effect resulting from 
A3/2 < Ai should in principle be observable in present 
experiments. 

To this end, let us compare the h = pair-breaking 
energies in an even and an odd grain, defined as the 
energy per electron needed to break a single pair at h = 
by flipping a single spin: for an even grain, it is Q e = 
7}{Ei — E ) h=0 , i.e. simply half the spectral gap discussed 
above; for an odd grain, it is Q = ^(£3/2 — Ei/ 2 )h=o- 

Within mean-field BCS theory, one would evaluate 
these using the same pairing parameter A for all states 
and [(ej — fi p ) 2 + A 2 ] 1 / 2 for the quasi-particle excita- 
tion energies associated with having the single-particle 
state \j, ±) definitely occupied, or empty, with parity- 
dependent chemical potential 2 ] fi p — eo — pd/2. This 
would give flf cs = [(d/2) 2 + A 2 } 1 / 2 an d qBCS = 
[d 2 + A 2 } 1 / 2 , implying that the difference flf cs - nf cs 
is strictly > (Fig. ||, dotted lines). For d/A — > 00 this 
difference reduces to d, which is simply the difference in 
the kinetic energy cost required to flip a single spin when 
turning |i) into ||)o (namely 2d), relative to that when 
turning |0)o into |l)o (namely d). 

In contrast, using the present theory to go beyond 
mean-field BCS theory, one finds numerically that £l > 
Sl e only for sufficiently large level spacings (d/A > 0.6, 
see Fig. ^, dotted lines); for smaller d one has il < fl e , 
implying that it costs less energy to break a pair in an odd 
grain than an even grain, even though the kinetic-energy 
cost is larger (2d vs. d). This happens since A 3 / 2 < Ai, 
which reflects a parity effect caused by pair-blocking by 
the extra unpaired electron in |3/2) relative to |1). The 
theoretical result that fl /Q e < 1 for sufficiently small 
d/A can be viewed as a "pair-breaking energy parity ef- 
fect" which is analogous to the "ground state parity 
effect" A1/2 < Ao, but which, in contrast to the lat- 
ter, should be observable in the experimentally available 
fixed-N eigenspectra. 

What are fl e and f2 Q in RBT's experiments? Unfor- 
tunately, the present data do not give an unambiguous 
answer: on the one hand, the h = data allow the deter- 
mination of f2 e = 0.25meV [half the h = energy differ- 
ence between the two lowest lines of Fig. ||(a)], but not 
of fl , since breaking a pair is not the lowest-lying exci- 
tation of an odd system at h = (which is why Fig. ^|(b) 
has no spectral gap). On the other hand, both f2 e and 
Q can be found from / data, since by Eq. ( |l9| ) 
they are equal to the level-crossing fields ft.0,1 = and 
^1/2,3/2 = ^o, whose values were deduced from the ex- 



f2 e -value somewhat smaller than the above-mentioned 
0.25meV determined at ft. = 0. The reasons for this dif- 
ference are presumably (i) that the actual ^-factors are 
not precisely 2 (as assumed), and (ii) that the experi- 
mental spectral lines are not perfectly linear in h (having 
a small /i 2 -contribution due to orbital diamagnetism, ne- 
glected in our model). 

Nevertheless, if we assume that these two complica- 
tion will not significantly affect the ratio ft.1/2, 3/2/^0,1 
(since ft-1/2,3/2 an d ^0,1 presumably are influenced by 
similar amounts), we may use it to estimate the ratio 
fl /fl e = 4.25/4 = 1.06 ± 0.1. This ratio is slightly 
smaller than that expected from the mean-field BCS ra- 
tio QBCS/qBCS ^ n at d j]± ^ .63 5 i.e. consistent with 

the pair-breaking energy parity effect. However, the dif- 
ference between 1.06 and 1.1 is probably too small to 
regard this effect as having been conclusively observed. 

We suggest that it should be possible to conclusively 
observe the pair-breaking energy parity effect in a some- 
what larger grain with /11/2, 3/2 < ^0,1 (implying f2 /f2 e < 
1), i.e. in Regime II of Fig. ||. (This suggestion assumes 
that in regime II the complicating effect of orbital dia- 
magnetism is still non-dominant, despite its increase with 
grain size.) To look for this effect experimentally would 
thus require good control of the ratio d/A, i.e. grain size. 
We suggest that this might be achievable if a recently- 
reported new fabrication method, which allows system- 
atic controL of grain sizes by using colloidal chemistry 
techniquest-3, could be applied to AI grains. 



3. Parity Effect in the Limit d/A 3> 1 

Since the parity effects discussed above are based on 
the observation that the amount of pairing correlations, 
as measured by A s , have a significant s-dependence, they 
by definition vanish for d > d c $, because then A s = 
for all s. Matveev and Larkin (ML)La have pointed out, 
however, that there is a kind of parity effect that persists 
even in the limit d/A 3> 1, which in the present theory 
we would call the "uncorrelated regime" (since there the 
A' defined in Eq. (Q) would be <§; A): when one extra 
electron is added to an even grain, it does not participate 
at all in the pairing interaction, simply because this acts 
only between pairs; but when another electron is added 
so that now an extra pair is present relative to the initial 
even state, it does feel the pairing interaction and makes 
a self- interaction contribution —Xd to the ground state 
energy. To characterize this effect, they introduced the 
pairing parameter 



a ML 



E N+1 



E^ 



2 ). 



(23) 



pcrimcntal tunneling spectra in Section [V B . This yields 
tt e = 0.23 ± O.OlmeV and n o = 0.24 ± O.OlmeV, i.e. a 



with A^=even. In first order perturbation theory in A, i.e. 
using Ep^ P = o(p\H\p)o (where |p)o is the uncorrelated 
Fermi ground state with N + p electrons) , one obtains 
A„ :P ° rt = ^Xd. This illustrates that this parity effect 
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exists even in the complete absence of correlations, and 
increases with d. 

Since our variational ground states \p) reduce to the 
uncorrelated Fermi states \p)o when Ap = 0, the above 
perturbative result for d/A can of course also be retrieved 
from our variational approach: we approximate Eq and 
-^1/2 1 k v £o(d) and £ l / 2 (d), respectively, both of which 

were calculated above, and E^ +2 by £o(d) — Xd, since it 
differs from £o{d) only by an extra electron pair at the 
band's bottom, whose interaction contribution in Eq. (|^) 

is — Ad(uj s ^) 4 . Thus the variational result for ML's parity 
parameter is 



A 



£ 1/2 (d) - £ (d) + Xd/2 , 



(24) 



(see Fig. [10|), which reduces to the perturbative result 
Ap L,pelt for d > d c ,o- The reason why this parity effect 
did not surface in the discussions of previous sections in 
spite of its linear increase with d is simply that there 
we were interested in correlation energies of the form 
£ — £° in which effects associated with "uncorrelated" 
states were subtracted out [see e.g. Figs. |(b) and |(b)] . 

The perturbative result Ap L ' pert = i Ad is in a sense 
trivial. However, ML showed that a more careful calcula- 
tion in the regime d/A » 1 leads to a non-trivial upward 
renormalization A of the bare interaction constant A, 



A = 



A 



1 - Xlog(ujc/d) 



(25) 



To obtain Ap L with logarithmic accuracy, A in A^L.pert 
is replaced by this renormalized A, with the result 



A ML 



dj (2 log d/A). 



(26) 



(The range of validity of this result lies beyond that 
shown in Fig. [l(], which therefore does not also show 
Ap L .) This logarithmic renormalization, which is be- 
yond the reach of our variational method (but was con- 
firmed using exact diagonalization in Ref. p7f ), can be re- 
garded as the "first signs of pairing correlations" in what 
we in this paper have called the "uncorrelated regime" 
[in particular since | Ap L | increases upon renormalization 
only if the interaction is attractive, whereas it decreases 
for a repulsive interaction, see Eq. (p5|)]. 

Unfortunately, Ap L is at present not measurable, for 

the s ame ex perimental reasons as apply to Eq°, see Sec- 
tion tvcii. 



D. Time Reversal Symmetry 

When defining our model in Eq. (^), we adopted a 
reduced BCS Hamiltonian, in analogy to that conven- 
tionally used for macroscopic systems. In doing so, we 
neglected interaction terms of the form 



d K i ^,i',3 , )4,+ C ]- C i'- C . 

iji'j' 



(27) 



between non-time-reversed pairs c| + cj_ , following An- 
derson's argumeniS that for a short-ranged interac- 
tion, the matrix elements involving time-reversed states 
Cj + c^_ are much larger than all others, since their or- 
bital wave-functions interfere constructively^. Interest- 
ingly, the experimental results provide strikingly direct 
support for the correctness of neglecting interactions be- 
tween non-time-reversed pairs of the from (|27j ) at h = 0: 
Suppose the opposite, namely that the matrix elements 
X(j + k, j, j' + fc, j') were all roughly equal to A for a finite 
range of k- values (instead of being negligible for k ^ 0, as 
assumed in H te &), Then for 2s < k, one could construct 
a spin-s state \s)' with manifestly lower energy {£') than 
that {£) of the state \s) of Eq. (|l5[): 



i s )' = n 4 + n o 



,(s) 



0) t 
■>.- c 



(i+2s) 



+ 4_)|Vac). (28) 



3 = -m 



Whereas in \s) pair- mixing occurs only between time- 
reversed partners, in \s)' we have allowed pair- mixing be- 
tween nora-time-reversed partners, while choosing the 2s 
unpaired spin-up electrons that occupy their levels with 
unit amplitude to sit at the band's bottom (see Fig. pi] ). 
To see that |s)' has lower energy than |s), 



£'s 



£' 



£'" < £ c a 



(29) 



we argue as follows: Firstly, £'° s = since the corre- 
sponding uncorrelated states |s)g and |s)o are identical 
[and given by Eq. @]. Secondly, A' s = A (> A s ), and 
hence £' c ° rr = £^°^\< £f rr < 0), because the 2s un- 
paired electrons in |s)' sit at the band's bottom, i.e. so 
far away from ef that their blocking effect is negligible 
(whereas the 2s unpaired electrons in |s) sit around Ef 
and cause significant blocking). Thus Eq. ( p9| ) holds, 
implying that |s') would be a better variational ground 
state for the interaction (p7| ) than |s). 

Now, the fact that £' c s = £^ orr is independent of s 
means that flipping spins in |s)' does not cost correlation 
energy. Thus, the energy cost for turning |0)' into |1)' 
by flipping one spin is simply the kinetic energy cost d, 
implying a threshold field h' j = d/2 [see Eq. (|i~9|)]; in 
contrast, the cost for turning |0) into |1), namely 2f2 e , 
implies a threshold field /io,i — Q e , which (in the regime 
d~A) is rather larger than d/2. The fact that RBT's 
experiments [Fig. 0(b)] clearly show a threshold field /io,i 
significantly larger than d/2 shows that the actual spin- 
1 ground state chosen by nature is better approximated 
by |1) than by |1)', in spite of the fact that £' \ < £\. 
Thus the premise of the argument was wrong, and we 
can conclude that those terms in Eq. (|27|) not contained 
in H re d can indeed be neglected, as done in the bulk of 
this paper. 
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V. CONCLUSIONS 

Citing the extensive literature in nuclear physics on 
fixed- N projections of BCS theory, we argued that a rea- 
sonable description of ultrasmall grains is possible using 
grand-canonical BCS-theory, despite the fact that such 
grains would strictly speaking require a canonical de- 
scription. Using a generalized variational approach to 
calculate various eigenenergies of the grain, we demon- 
strated the importance of the blocking effect (the re- 
duction of pair-mixing correlations by unpaired elec- 
trons) and showed that it becomes stronger with de- 
creasing grain size. The blocking effect is revealed in the 
magnetic-field dependence of the tunneling spectra of ul- 
trasmall grains, in which pairing correlations can be suf- 
ficiently weak that they are destroyed by flipping a single 
spin (implying "minimal superconductivity"). Our the- 
ory qualitatively reproduces the behavior of the tunneling 
thresholds of the spectra measured by Ralph, Black and 
Tinkham as a function of magnetic field. In particular, it 
explains why the first order transition from a supercon- 
ducting to a paramagnetic ground state seen in thin films 
in a parallel field is softened by decreasing grain size. Fi- 
nally, we argued that a pair-breaking energy parity effect 
(that is analogous to the presently unobservable ground 
state energy parity effect discussed previously) should be 
observable in experiments of the present kind, provided 
the grain size can be better controlled than in RBT's 
experiments. 
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APPENDIX A: ANALYTICAL LIMITS 
1. d — > and Euler-MacLaurin Expansion 

When the level spacing d tends to zero the theory re- 
duces to the conventional BCS variational and mean field 
approach. We can calculate the properties of a supercon- 
ducting system to first order in d by expanding the BCS 
solution around d = 0. In doing so, we focus on the 
ground states \s) of each spin-s sector of Hilbert space. 

While in the bulk limit (d = 0) the shift -Xd(v^) 2 
in the single-electron energies £j just after Eq. (O) is 
unimportant, it influences the behaviour of an ultrasmall 
grain by effectively increasing the level-spacing near the 
Fermi surface. Its effect is largest for s = 0, since for 



s ^ the states at the Fermi surface, where the devi- 
ation of v^' a%> from or 1 is largest, are blocked. For 

simplicity we neglect the ^'"''-dependence in £j in the 
following calculation, using £j = Ej — fi — \d0(—(ej — fi)), 
and therefore good agreement with numerics can only be 
expected for d <§C A and s / 0. Within this approxi- 
mation for £j, fi lies halfway between the topmost dou- 
ble occupied and lowest completely empty level in \s)o- 
fj, = EQ — d(6 P! o + \)/2. Note that /i does not lie exactly on 
one of the levels in the odd case (p — 1) as one might have 
expected at first sight, but halfway between the topmost 
doubly-occupied and lowest completely empty level. 

We shall calculate the pairing parameter A s (d) in the 
small-d limit by calculating the first terms of its Taylor 
series: 



(Al) 



A s (d)-(l + dd d + ^-d 2 d )A s (0) 



To this end, it suffices to solve the gap equation rtl2), as 
well its first and second derivatives with respect to d, for 
d = 0. This can be done by rewriting Eq. ( |l2| ) using the 
Euler-MacLaurin summation formula, 



i/\=dj2f(jd)~ r de/co 

i— nr. J jtid 



(A2) 



+^[f(jod) + f(nd)} + ^[f'(jod) + find)}, 

with f(jd) = [(jd) 2 + A*]" 1 /* jo = s + (1 + A)/2 and 
ji = uj c /d. The s-dependence has now been absorbed in 
the lower bound jo of the sum. The negative branch of 
the sum is identical to the positive since fi lies halfway 
between the topmost doubly-occupied and lowest com- 
pletely empty level. It therefore suffices to calculate the 
positive branch times two. Setting d = in Eq. (A2) 
yields the well-known BCS bulk gap equation, whose so- 
lution is, by definition, A s (0) = A. The first and second 
total d-derivatives of Eq. (A2) yield ddA s (d = 0) = — s 
and djA s (d = 0) 
from Eq. ( |Al|) is 



-s /A, so that the desired result 



A s (d) ~ A - (s + X/2)d - 



(s + X/2) 2 d 2 
2A~ 



(A3) 



We next calculate the eigenenergies £ s by evaluating 
Eq. (||) up to first order in d, where the sums again are 
evaluated with the help the the Euler-MacLaurin for- 
mula. Since we are interested in the effects of pairing 
correlations we subtract the energy £® of the uncorre- 
cted Fermi sea \p)o- 



(even) £ s - 
(odd) £ s -£l 



-t r + (l + f)AA + 2s A- 
" 2 i + ^Xsd; 



Al 

2d 

-Or 
_Ai 

2d 

-(s 2 



fAA 
4 



4 

2sA 
^ As 



(A4) 



^Xs + |)d. 
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The d 1 term is the bulk correlation energy, which is 
slightly renormalised by the intensive (1 + ^AA-term, 

which in turn stems from the v 4 -terms of Eq. (JsJ) . 2s A is 
the bulk excitation energy for 2s quasi-particles. The 
d -term is the first-order correction for discrete level- 
spacing. 



2. d near d c and the Small Delta Expansion 

The other analytically tractable limit is d ^ A s , which 
holds for d near the critical spacing d c s where A s van- 
ishes. 

First, we derive an expression for the critical d c . s by 
solving the gap equation with vanishing pairing parame- 
ter A s for d: 

T = E - = *(<"c/dc. + 1) - *0'o)- (A5) 



\l/(x) denotes the Digamma function and jo equals s - 
ijp again. Remembering that A = l/log(2^) an 

exp(5'(x)) ~ x — \ for large x this equation reduces to 



log ('H^i ) = -*(s + i±A) (A6) 



A 

4, s = ^exp(-*( S +i±A)). 

For s > 1 this can be simplified to 

A 



2s + A 



(A7) 



(A8) 



Numerical values of d c<s J A (A = 0.224) are 
2.36,0.77,0.44,0.31,.. . for s = 0, ±, 1, f, . . . respectively. 
Near d CjS the pairing parameter vanishes like 



A s ~ A J 1 - - — for d > A s and s > 0, (A9) 

\/ Clr- c 



which we shall now show. 

Since for the spin-s ground states with vanishing pair- 
ing parameter electron and hole pairs are symmetrically 
distributed around the Fermi surface, Eq. (13) again 
yields lu = £o — d(S Pt o + A/2). We turn to the gap equa- 
tion ([l2|). The spin dependence has been absorbed in j . 
The positive and negative branches of the restricted sum 
are identical (because of the special symmetric value of 
fx), with |£| ranging from d(s + i±^) = dj to lu c . It 
therefore suffices to calculate the positive branch times 
two: 



1 



uj c /d 
3= Jo 



Al/d* 



-1/2 



io c /d c s 

E 

j= s +(l + A)/2 



uj c /d 

E 

j= s +(l+A)/2 
u) c /d 

E 

j= s +(l+A)/2 



A 



2d?f 

A 2 
2d 2 f 



(A10) 



To obtain Eq. (A10), the square root was expanded using 
A s <C d. The remaining sums can be expressed by the 
polygamma functions using the identity 



n 1 



fe=i 



(m - 1)! 



(All) 



Replacing the sums by the Polygamma functions and col- 
lecting terms leads to 



A|_r 

Ad? 



Now assume that d is close to d r 



1) = 
1 + A 



(A12) 



Sd and 



(5d <C d c s . Expand the left hand side in Sd and use the 
asymptotics for iff' (on the left hand side) and ^" (on the 
right hand side) for the large oj c /d argument. Also the 



*"(s 



-term is approximated by its asymptotic form 



Sd 



Ad 2 



*"(s + ±) 



A 2 = 4d 2 s 2 ^4- 



(A13) 
(A14) 

A s =AJl--^-. (A15) 
The last step was performed by remembering that 



Ad 2 



Ad 2 



A 2 for s ^ 0. 



Although Eq. (A9) was derived for d near d CjS , it turns 
out to have a surprisingly large range of validity: its 
small-d expansion in powe rs o f d/A agrees (at least) up 
to second order with Eq. (A3), and for s > 1 it in fact 
excellently reproduces the numerical results for A s (d) for 
all d. 

For s = the asymptot ic ex pansion of 'J'" breaks down. 
Therefore directly from ( A13| ) we deduce 



A f 



Ad 2 d r 



12.1 d r 



(A16) 



where we used $"(^4= 



— 12.1. This result gives good 
agreement with numerics near d c . s= Q, but obviously has 
the wrong d — > limit. 
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APPENDIX B: I-V CHARACTERISTICS OF AN 
ULTRASMALL NSN SET 

The I-V characteristics of a SET with an ultrasmall 
superconducting grain as island, i.e. an ultraspall NSN 
SET, were examined by Tichy and von Delftnj. They 
described the discrete pair-correlated eigenstates of the 
grain using the parity-projected mean-field BCS theory of 
Rcf. [l^. Although this approach is too crude to correctly 
treat pairing correlations of excited states (since for all 
even (or odd) ones the same Ao (or A1/2) is used), it 
does treat the even and odd ground states correctly. It 
therefore enables one to understand how the odd-even 



ground state energy difference Eq £ = (E 1 / 2 
influence the SET's I-V characteristics. 

Using tunneling rates given by Fermi's golden rule and 
solving an appropriate master equation, Tichy calculated 
the tunnel current through the SET as a function of 
transport voltage V and gate voltage V g at zero magnetic 
field. In an ideal sample, the I-V characteristics are 2e- 
periodic in the gate charge Q g = V g C g + const; one such 
period is shown in Fig. n3. The usual Coulomb-blockade 
"humps" centered roughly around the degeneracy points 
Q g /e — 2m±i are decorated by discrete steps, due to the 
grain's discrete eigenspectrum. In RBT's experiments V g 
was fixed near a degeneracy point and the current mea- 
sured as function of V (for a set of different H- values). 
When following a line parallel to the U-axis in Fig. [12, 
the positions of the steps in the current thus correspond 
to the H = eigenenergies of RBT's tunneling spectra 
in Fig. |. 

The reason for 2e- instead of e-periodicity are pairing 
correlations: Firstly, the grain's odd-even ground state 
energy difference E°q & causes a shift in the degeneracy- 
point values for Q g /e from 2m±i to 2m±(| + Eq ' e / Ec) ■ 
Secondly, tunneling spectra measured in the U-direction 
in Fig. [lj show a plateau after the first step if the fi- 
nal state after tunneling is even (i.e. for i + E° G /e /E c < 

Qg/e < § — Eq/Ec), but not if it is odd, correspond- 
ing to the presence or absence of a large spectral gap in 
the tunneling spectra of Fig. ^(a) or (b); this is due to 
the energy cost to break a pair, and the plateau's width 
is simpl y twic e the even pair-breaking energy 20 e (see 
Section jlVC 21). 
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FIG. 1. A cartoon depiction of why "superconductivity 
breaks down" when the sample becomes sufficiently small. 
Each vertical line represents a pair of single-particle state | j±) 
with energy £-,-, for three different mean level spacings d, cor- 
responding to (a) a "large" grain (d <C A); (b) a "small" 
grain (d ~ 0.25 A); (c) an "ultrasmall" grain (d ~ A). In all 
three plots, the height of each vertical line equals the func- 
of standard bulk BCS theory, illustrat- 



tion u 2 v 2 



l A 

4 e' 2 + A 2 



ing the energy-regime (of range A around sf) within which 
electrons are affected by pairing correlations. Loosely speak- 
ing, the number of single-electron states A/d in this regime 
corresponds to "the number of Cooper pairs" of the system. 
Evidently, when d/A ~ 1 as in (c), "the number of Cooper 
pairs" becomes less than one and it no longer makes sense to 
call the system "superconducting" . 
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FIG. 3. Properties of spin-s ground states |s) [compare 
Eq. (p^)] : (a) The pairing parameters A s (d)/A for some 
spin-s ground states \s), as a functions of d/A. The criti- 
cal level spacings d CjS at which A a (d c , s ) = are found to be 
2.36, 0.77, 0.44, 0.31, ... for s = 0, 1/2, 1, 3/2, . . ., respectively, 
(b) The energy densities (f s — f^ 2 )d/A 2 (solid lines), plotted 
as functions of d/A for h = 0, of some pair-correlated spin-s 
ground states \s) relative to the uncorrelated spin-p/2 Fermi 
sea |p/2)o, and for comparision the relative energy densities 
(£° — £p/ 2 )d/A 2 (dashed lines) of the corresponding uncorre- 
lated paramagnetic states |s)o (obtained from |s) by setting 
A s = 0) . We call the plotted quantities energy densities since 
the normalization factor d/A 2 contains d ~ Vol -1 . The solid 
and dashed spin-s lines meet at the critical level spacing d c<a , 
above which no pairing correlations survive (so that the rela- 
tive energy densities equal (s 2 —p/4+(s—p/2)\)d 2 / A 2 there). 



FIG. 2. Two examples of states in the spin-| sector of 
Hilbert space: (a) the ground state ||) and (b) the excited 
state 1 1 , 2) . The single-particle levels are drawn at h = 0, and 
we indicated schematically how states are paired according to 
(m +Vic\ + c\_) in the BCS- like Ansatze (|l|) and (jl|) for |§) 
and 1 1 , 2) , with solid or dashed ellipses connecting states that 
would be completely filled or empty in the absence of pairing 
correlations. 
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1 2 

Level spacing d/A 



FIG. 4. Properties of excited spin-i states \ [compare 
Eq. (p7|)]: (a) The pairing parameter for some spin-i 

states (j = 0, ...,4), together with Ao of the spin-0 

ground state |0) (the outermost curve). The larger j, the 
closer A-i/2,j approaches the spin-0 value Ao. (b) The relative 
energy densities (£1/2,3 ~ ^1/2, o)d/A 2 (solid lines) of 
relative to ||,0}o = |^}o, and for comparison the relative 
energy densities (£i/ 2 j ~ £1/2 o)^/^ 2 (dashed lines) of the 
corresponding uncorrelated state j)o- For excited states 
the solid and dashed lines meet at a larger d than for the 
ground state, i.e. in excited states pairing correlations survive 
down to smaller grain sizes than in the corresponding ground 
state. 
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FIG. 6. Experimental tunneling spectra measured by RBT 
(Fig. 3 of Ref. m. The distances between lines give the fixed- N 
excitation spectra of (a) a n even and (b) an odd grain, as 
explained in Section IV A. The vertical dashed lines indi- 
cate the first four level-crossing fie lds H s , s i (assigned by com- 
parison with Fig. \A see Section 



tVB), namely #0,1 = 4T, 



Hi 



/2,3/2 : 



= 4.25T, #i,2 = 5.25T and H- 



3/2,5/2 



= 6.5T with un- 



certainty ±0.13T (half the Jf-resolution of 0.25T). 




FIG. 5. The level-crossing fields h B>s ,(d)/A [see Eq. @] 
for the cascade of ground state changes (GSCs) (so,si); 
(si,S2); ... that occurs as h increases from at given d. 
Some lines are labeled by the associated GSC (s, s') (where 
£ s i drops below £ B as h increases past /i SjS '). (Level crossing 
fields not associated with a GSC are not shown.) The order 
in which GSCs can occur within a cascade (i.e. the order of 
h SyS i lines encountered when moving vertically upward in the 
figure) depends sensitively on d, and an infinite number of 
distinct regimes (cascades) I, II, III, . . . can be distinguished. 



19 



a) 



> 
en 

CD 
LU 

b) 




0.5 

Magnetic Field h/A 



FIG. 7. The theoretical odd-to-even and even-to-odd tun- 
neling spectra (AEif — AE m in{Q)) / A predicted for an ul- 
trasmall superconducting grain as a function of magnetic 
field h, for two different level spacings: (a) d = 0.67A 
and (b) d — 0.34A (corresponding to regimes I and III of 
Fig. 0, respectively). Some lines are labeled by the cor- 
responding Si — » s'i tunneling transition. Not all possible 
higher lines (corresponding to excited final states \s,j)) are 
shown. Vertical dashed lines indicate those level-crossing 
fields h s iS i [see Eq. ( |l9| at which kinks or jumps oc- 
cur, with h ,i < fti/2,3/2 < hi t 2 < ^3/2,5/2 in (a) and 
hi/2, 3/2 < h , 2 < h 2 ,3 in (b). 
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Section 



9. Parity effect for the 
h (Ei - E )h=o and £l a 



pair-breaking energies 
1(^3/2 - E 1/2 ) h =o (see 
IV CS): when calculated naively using conventional 
mean-field theory (dashed lines), the pair-breaking energies 
obey f2 > f2 e for all d/A; in contrast, when calculated within 
generalized variational BCS theory (solid lines), Q < f2 c 
for d/A < 0.6; this reflects a parity effect, namely that 
A3/2 < Ai, which is caused by the extra unpaired electron in 
1 3/2) relative to |I). 
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FIG. 10. r^he parity parameter Ap L discussed by Matveev 
and LarkinEJ, calculated perturbatively for the uncorrelated 
Fermi sea (/\^ L >P crt — I / \ c ( i dashed line), and using our gen- 
eralized variational BCS-approach (Ap L,var of Eq. (|2^), solid 
line) . The renormalized result Ap L of Eq. ( pj|) given by ML 
is not shown since its range of validity is d/ A 3> 1. 



FIG. 8. The first jump in the lowest line of a tunneling 
spectrum occurs at the level-crossing field hcc(p,d) = h SOt31 
that induces the first ground state change (so,si). The size 
of this jump in units of A, namely |Ai5 sl j/ — Ai? SOi /|/A, 
is shown here as function of d/A for e — > o (solid line) and 
o — *• e (dashed line) tunneling spectra. Both lines approach 
the CC value 1-1/ V% = 0.29 as d -> 0. Their non-monotonic 
behavior is due to the discreteness of the level spacing; the 
kinks occur at the regime boundaries of Fig. [| indicated here 
by vertical dashed lines. 
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FIG. 11. Schematic represen- 

tations of the non-time-reversed-pairing state 1 3/2) ' defined 
in Eq. ( pj| ) . The energies Ej =F h of the single-particle states 
\j, ±) are drawn (a) for h — and (b) for 2h = 3d. We indi- 
cated schematically how non-time-reversed states are paired 
according to (m + u i c ( i+3 ) + c l_) m t ne BCS-like Ansatz (pi[), 
with solid or dashed ellipses encircling states that would be 
completely filled or empty in the absence of pairing correla- 
tions. 




FIG. 12. I-V characteristics for a SET with an ultrasmall 
superconducting grain as island (from Ref. |3rj| ) . The current is 
plotted as a function of gate charge (Qg/e) and transport volt- 
age (eV/Ec)- Pairing correlations shift the degeneracy-point 
values of Q g /e away from their e-periodic values of 2m ± i 

by ±Eq £ / Ec (see text). To better reveal the figure's charac- 
teristic features, it was plotted using a ratio Eq / Ec — 0.1, 
very much larger than the typical values of < 0.01. 
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